Collisional de-excitation in a quasi-2D degenerate Bose gas 
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We separate a Bose-Einstein condensate into an array of 2D sheets using a ID optical lattice, and 
then excite quantized vibrational motion in the direction normal to the sheets. Collisions between 
atoms induce vibrational de-excitation, transferring the large excitation energy into back-to-back 
outgoing atoms, imaged as rings in the 2D plane. The ring diameters correspond to vibrational en- 
ergy level differences, and edge-on imaging allows identification of the final vibrational states. Time 
dependence of these data provides a nearly complete characterization of the decay process including 
the energies, populations, and lifetimes of the lowest two excited vibrational levels. The measured 
decay rates represent a suppression of collisional de-excitation due to the reduced dimensionality, a 
matter wave analog to inhibited spontaneous emission. 

PACS numbers: 03.75.Kk, 05.30.Jp 
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Most quasi-2D quantum systems have been realized 
with electrons in semiconductors, where a ID potential 
confines the electrons to the lowest quantized vibrational 
states in one direction, while leaving them essentially free 
in the remaining two. Recently it has become possible to 
confine degenerate atomic Bose gases to 2D and 
investigate vibrational excitations in the tightly confined 
direction Trapped 2D atomic gases provide experi- 
mental opportunities unavailable in electron systems. For 
example, unlike semiconductors, the atomic system is 
nearly defect-free. Further, the dynamic control of the 
confining potential, coupled with an ability to image the 
atoms, enables the direct detection of the excited state 
population and the momentum distribution. Quantized 
vibrational states are an ingredient in proposals to real- 
ize exotic states of matter, such as striped or super-solid 
phases 0, 0, 0| , and are possible motional qubit states 
for quantum computation Stronger confinement 

(beyond that described herein) can also change the na- 
ture of collisions [Hill [12. 

Here we study the vibrational relaxation of a quasi- 
2D degenerate Bose gas where quantized motion in the 
tightly confined direction plays a role analogous to an 
internal degree of freedom. We transfer a large fraction 
of atoms into excited vibrational states, creating highly 
non-equilibrium atom populations. In this system, atom- 
atom collisions provide the only significant relaxation 
mechanism, transferring "internal" energy to 2D kinetic 
energy. We directly observe atom populations as outgo- 
ing rings (Fig. ^) representing distinct decay channels. 
The excited state lifetimes are enhanced due to the re- 
duction in the density of final scattering states relative 
to scattering in an unconfined 3D gas. 

We create independent 2D sheets (or pancakes) of 
atoms by applying a deep ID optical lattice to a 3D 
Bose-Einstein condensate (BEC). In the tight direction 
(z) the system is well-described by a single-particle ID 
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FIG. 1: Absorption images of vibrationally excited atoms af- 
ter a 1 ms decay and subsequent TOF, with 60% initial pop- 
ulation in the a) n = 1 and b) n = 2 vibrational levels. In 
the top views, outgoing rings (txcF ~ 7-1 ms) correspond to 
different in-plane energies imparted to the atoms from various 
decay channels. Viewed from the side (tTOP ~ 13.1 ms), the 
rings appear as rectangles, where the vertical momentum dis- 
tribution identifies the final vibrational states. The schematic 
of the side view shows the expected distributions from the 
dominant decay processes to n = 1 (white rectangles) and 
n = (hatched rectangles). 



Schrodinger equation, yielding discrete vibrational levels 
labeled by an index n. Atoms, Raman-excited from the 



2 



ground state n = 0ton = lor2, collide and decay. By 
imaging after a time-of-flight ixOF we identify the mo- 
mentum and population of atoms in the various final vi- 
brational states. We extract excited state lifetimes from 
time sequences of these single-shot vibrational spectra. 

We produce a magnetically trapped ^''Rb BEC with 
up to 2.5 X 10^ atoms in the \F = l^mp = —1) state [is). 
The BEC is separated into a stack of about 80 pancakes 
by an optical lattice with period d = 410(1) nm 0| ver- 
tically aligned along z When in the ground vibra- 
tional state, the largest pancake has N « 4.6(10) x 10'^ 
atoms with a 2D chemical potential /i2D ~ 1-6(2) kHz; 
the resulting peak 2D and 3D densities are 2.2(2) x 
10^ cm-2 and 2.9(3) x lO" cm'^, respectively. The 2D 
Thomas-Fermi radii are — 11.0(6) fim and Ry = 
12.0(7) /zm. In the combined magnetic and optical po- 
tential, the in-plane oscillation frequencies are lUx/2tt = 
55(1) Hz and ujy/2TT ~ 50(1) Hz. The lattice is raised 
in 200 ms with an exponentially increasing ramp (50 ms 
time constant). This timescale is chosen to be adi- 
abatic with respect to mean-field interactions and vi- 
brational excitations. By pulsing the lattice and ob- 
serving the resulting atom diffraction [l^ . we measure 
a lattice depth of s = 77(4), expressed in units of 
Er = K^/Snuf = h X 3.42(2) kHz. For a single well 
of the deep sinusoidal potential, the energy spacings are 
En+i—En — Er[2^— [n + l)] (these energies, which in- 
clude the lowest order anharmonic correction, are shifted 
at the 1% level by the inter-atomic interaction). The 
harmonic frequency is uiz/Itt = 2Euy/s/h = 60(1) kHz. 

Raman transitions between vibrational levels [l^ are 
driven by a pair of laser beams. The nearly counter- 
propagating Raman beams are oriented approximately 
along z [13, are detuned 82 GHz below the *'^Rb D2 
transition, and have a relative detuning d ranging from 
50 to 120 kHz. A 1 ms Raman pulse excites a fraction 
/„ of the atom population to either n = 1 or n = 2. 

The anharmonicity of the potential allows us to selec- 
tively transfer populations between desired vibrational 
levels, provided the pulse duration tp is long enough 
that its Fourier spread resolves the > Er/H difference 
from unwanted transitions, tp must also be shorter than 
the vibrational lifetime r. For our experiment (where 
T }t }t fi/Eii), we find that by detuning below Ra- 
man resonance and chirping 15 kHz through resonance 
in 1 ms, we controllably transfer up to 65% of the atoms 
to either n = 1 or n = 2 (simulations indicate a maxi- 
mum transfer of around 75% in this situation) without 
significantly populating unwanted states (~ 5%). 

The vibrationally excited atoms are allowed to decay 
for variable hold times thoid ranging from 10 /iS to 10 ms; 
then the lattice is turned off in 200 /iS. The magnetic 
trap is then turned off in « 300 fis, and the atom cloud 
expands for ^tof before absorption imaging. The re- 
sulting atomic distribution, which constitutes an average 
over the ~ 80 pancakes, is imaged along one of two axes: 



Ring 


Decay path 


Ef — Ei 


rab 


Central ring 


11 - 


^ 20 


Er/2 


0.0243 




22 - 


^ 31 


Er/2 


0.0257 




22 - 


^ 40 


2Er 


0.0007 


Inner ring 


11 - 


^ 00 


huj^ - Er 


0.2202 




22 - 


11 


huoz — 2Er 


0.1556 




22 - 


^ 20 


hu^ - ZErI2 


0.0020 




20 - 


^ 00 


hw^ - 3Er/2 


0.1128 


Outer ring 


22 - 


^ 00 


2huj^ - 3Er 


0.1075 



TABLE L Release energies (including the lowest order anhar- 
monic corrections) and matrix elements for the experimen- 
tally relevant decay paths and lattice depth s = 77 (sorted in 
order of increasing decay energy) . The energies denote the per 
atom difference between initial and final vibrational energies. 



perpendicular to the pancakes (top view, along z) and 
edge-on (side view, along y). 

The decay of a pair of atoms from initial vibrational 
states a and b to final states c and d (denoted ab cd) 
can be viewed as a 2D two-body inelastic collision pro- 
cess, where vibrational states in the confined direction 
play the role of internal states. The release of vibrational 
energy leads to back-to-back atom pairs with large mo- 
menta in the 2D plane. Figure H shows example TOF 
images from which we can infer the final momentum dis- 
tribution. Image coordinates are in units of recoil mo- 
mentum, hkji = 7r?i/d |2l| . and the corresponding recoil 
velocity is hkn/m = 5.6 mm/s. Each image pair (Figs. 
^ and^) contains nearly full spectral information of an 
excited state and its decay paths. 

The decay paths allowed by parity and conservation of 
energy are shown in Table For atoms excited to n = 1, 
only two decay paths contribute. The high-energy pro- 
cess, 11 — > 00, gives rise to the ring in Fig. QJi (top view), 
while the ring from the low-energy 11 ^ 20 process is not 
resolved from the central cloud. Each ring's radius corre- 
sponds to the atomic in-plane velocity and hence the per 
atom difference between initial and final vibrational en- 
ergies. In Fig. atoms excited to n = 2 decay through 
several channels. The processes 22 11, 20 ^ 00, and 
22 20 all contribute to the inner ring of Fig. QJd. Each 
process imparts an energy of about 2E}iy/s per atom; the 
differences in energy of Eii/2 are small compared to the 
ring's width and are unresolved in the image. The outer 
ring results solely from the 22 00 decay process. 

FiguresE^L and|2t) show the radial density of the atomic 
cloud after an angular integration of the data in Figs. ^ 
andQJ), respectively. We extract the vibrational energy 
spacings by fitting the data to Lorentzians. We attribute 
the atom background between the peaks to secondary 
scattering involving outgoing atoms and include it in the 
fits as a linearly sloping baseline. For atoms initially in 
n = 1 (Fig. [2Jl), the measured energy is /i x 55(2) kHz, 
in accord with the expected value for a 77(4) Er lattice: 



3 



Pi 




2 4 6 

Momentum [hk^] 



FIG. 2: Angular integrals of top view data, a) From Fig. 
the single ring gives rise to a peak located at fc = 4.0(1) ka. b) 
The two rings in Fig. ^ yield two peaks at = 3.9(1) kn and 
5.6(2) fefl. The diagrams schematically illustrate the decay 
channels which contribute to each peak. 



gives rise to atoms in the first BZ (hatched rectangle in 
Fig. ^), while 22 ^ 11 leads to atoms in the second BZ 
(white rectangles in Fig. ^). Note that the vertical cen- 
tral structure is taller since it contains atoms remaining 
in n = 2 (third BZ). 

We calculate short-time two-body branching ratios and 
decay rates in a single pancake using Fermi's golden 
rule. Due to tlie extreme anisotropy of our poten- 
tial the initial condensate wavefunctions can be approx- 
imated as a product of single particle wavefunctions, 
^'i(a;, 2:) « 'ip{x,y)ipni{z). Here it>(x, y) satisfies an ef- 
fective 2D Gross-Pitaevskii equation [lOj for n = atoms 
in the Thomas-Fermi limit; we assume ip{x,y) remains 
unchanged during the short duration of the experiment. 
iPm (z) solves the Schrodinger equation for the ID lattice 
potential and are nearly harmonic oscillator wavefunc- 
tions with an extent CTz = \/^lfnLOz- The final states 
^ f{x,y^z) = cxp[— ^(/cj;^ + fcy?;)](^„j: (z) are free particles 
in X and y which is justified since lOz ^ ^x,y- The rate 



for the scattering process ah - 
populations Na and Nh, and 



cd is Tf^NaNh, with atom 
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E^-Eii = ft. X 56.5(15) kHz. For the n = 2 data (Fig. |2b) 
we measure the energies of the outer and inner rings to 
be /i X 107(5) kHz and h x 52(2) kHz, respectively. The 
energies agree with the calculated values of E2 — Eq = 
h X 110(3) kHz and E2-Ei=hx 53.8(15) kHz {Ei is a 
weighted average using the matrix elements in Table [}. 
The branching ratio between the inner and outer rings 
estimated from Fig. |21 (ignoring the atom background) 
is 2.1(3). The value calculated from the matrix elements 
in TableQlis 2.9(3) for an excited fraction /2 — 0.60(5). 

The side view images in Figs. ^ and^ complement 
the top view images by identifying the final vibrational 
states. The 200 /iS lattice turn-off is adiabatic with re- 
spect to the 55 kHz vibrational frequency. As a result, 
the turn-off procedure maps quasi-momentum states in 
the ID lattice to the corresponding free particle momen- 
tum states .3|. For example, atoms in the n = 1 vi- 
brational state reside in the second Brillouin zone (BZ), 
and are mapped to a continuum of momentum states kz 
where kB.<\kz\ < 2kR ^22]. 

Figure^ shows atoms which were initially excited into 
n = \. In the side view, the atom cloud's extent along z 
reflects the mapping of quasi-momentum to momentum, 
and the extent in the horizontal direction reflects the final 
in-plane momentum. The ring in the top view appears 
as a rectangle in the side view (schematically illustrated 
by the hatched rectangle in Fig. ^). The dense, verti- 
cally aligned double-lobed structure at the center of the 
side view image is largely due to atoms which have not 
decayed and remain in the n = 1 state (second BZ). Fig. 
depicts atoms initially in n = 2. The process 22 00 



where /"J" = 2Tr \az J dzipa{z)(pb{z)ipc{z)(pd{z)\ is dimen- 
sionless (see Table Ht, and Ug = 5.3 nm is the s-wave scat- 
tering length 01 [23 ■ Parity considerations make = 
when a + b + c + d is odd. The total rate is a sum over 
energetically allowed final states c and d. 

For atoms excited to n = 1, the excited fraction /i(t) 
is governed by: 
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FIG. 3: Cross-sectional integral of Fig. ^1 (side-view). The 
vertical dashed lines delineate the lowest 3 BZs. For the 
specific example here, the fit (solid curve) gives fractions 
fo = 0.37, /i = 0.18, and /2 = 0.45 in each of the BZs. 
Inset; fhoid dependence of n = 2 population for a 60% initial 
fraction, showing the fit to the first 2 ms. 
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For atoms excited to n = 2, the fraction /2(<) is governed 
by the coupled equations: 



/2 



-27V[(r; 



-^22 



p20 

-•^ 00 



/o]/2 



/o = -27vr2;]/2/o 



These equations ignore small contributions from F^q and 
F40 and assume that once atoms decay, they leave the 
cloud and do not re-scatter. Parity considerations forbid 
the transition 10 — > 00 but not 20 00, so as /i 
the n = 1 decay rate per atom /1//1 0, but as /2 ^ 
the corresponding n = 2 rate is non-zero. 

The excited-state population as a function of ^hoid is 
found from a series of side-view images like those in Fig. 
n Integrating over x gives the atomic distribution along 
z; Fig. 13 is the z distribution corresponding to Fig. 
(ihoid = 1 ms). The dashed lines delineate the BZs. In 
this example, atoms are predominantly in n = and 
n = 2. To extract populations, we fit the data to a dis- 
tribution, flat within each BZ, convolved with a Gaussian 
(solid curve in Fig. The width of the Gaussian is fixed 
by applying this model to data in which only the ground 
state (first BZ) was occupied. Repeating the fitting pro- 
cess for different thoid yields the fractional population as 
a function of time (inset to Fig. |3J). We extract rates by 
fitting the first 2 ms (where we anticipate our model to 
be valid) to the expected solutions /i(t) and /2(i) H^l- 

The resulting rates per atom at t = are shown in Fig. 
01 For comparison we plot the predicted rates at t = 0, 
where the grey band reflects the experimental atom num- 
ber uncertainty. The agreement, with no free parameters, 
is good even though we neglected significant corrections 
from secondary collisions. The presence of atoms not on 
the rings (Fig. ^ and ^ indicates that secondary scat- 
tering is important. An analysis of Fig. [21 suggests that 
around 75% of outgoing atoms re-scatter, consistent with 
our theoretical estimates. These processes decrease the 
atomic density, and hence reduce the overall decay rate. 
(In separate experiments with lower density thermal sam- 
ples, we indeed observe significantly lower rates: below 
40 s~^.) The effect should be more significant for large 
initial excited fractions, a trend consistent with the data. 
The solid lines in Fig. 01 show the estimated effect of sec- 
ondary scattering from outgoing atoms which decayed 
during the 1 ms Raman pulse. 

In conclusion, we measured the collision-induced tran- 
sition rates between different vibrational states in a deep 
ID optical lattice using a new, single-shot, spectrographic 
technique. The reduced dimensionality of the final states 
suppresses the total scattering rate in a manner similar 
to the sup pression of spontaneous emission in a planar 
cavity |l7l |. In a separate calculation, we find that scat- 
tering into an unrestricted 3D space exceeds our observed 
rates by a factor between 2.5 and 3.5 (the variation stems 
from different possible assumptions in the calculation). 
Additional suppression of the rates would be achieved 
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FIG. 4: Initial decay rates, a) initial population in n = 1. b) 
initial population in n = 2. The thin lines represent the result 
of our model with no adjustable parameters, integrated over 
the distribution of pancakes expected from Thomas- Fermi ap- 
proximation. The grey regions reflect the variation in this 
model due to experimental number uncertainty. The thick 
solid lines include the estimated effect of secondary scatter- 
ing, using a 75% probability of a second event. 



by further modifying the final density of states lyL for 
example by confining the atoms into ID tubes |l3| . or 
in analogy with photonic band-gap materials, by apply- 
ing an additional in-plane lattice to open suitably placed 
band-gaps. The long lifetimes are expected to be useful 
in the context of producing correlated atomic systems. 
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